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ABSTRACT 

We analyzed the decays of B+ ^ D° + D+j{2S) and B+ ^ + D+j{lD) by naive fac- 
torization method and model dependent calculation based on the Bethe-Salpeter method, 
the branching ratios are Br{D^Djj{2S)) = (0.72 ± 0.12)% and Br{D^D^j{lD)) = (0.027 ± 
0.007)%. The branching ratio of decay B+ -f D^D+j{2S) D^D^K+ consist with the data 
of Belle Collaboration, so we conclude that the new state Dj^(2700) is the first excited state 
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BABAR Collaboration reported the observation of a new charmed meson called D^(2690) with a 
mass 2688 ± 4 ± 3 MeV and with a broad width T = 112 ± 7 ± 36 MeV [1], Belle Collaboration also 
reported a charmed state L»+j(2700), M = 2715 ± lltH MeV with width T = 115 ± 20^11 MeV [2], later 
modified to M = 2708 ± 9t\l MeV and T = 108 ± 23l|^ MeV [3]. Most authors believe that these two 
states should be one state, and the most possible candidate is the 2S or the S — D {2 ^Si — 1 ^Di) mixing 
cs 1- state [4, 5, 6, 7, 8, 9, 10, 11]. 

Recently, we have analyzed the decays of the Dfj{2S) and D^(ID), we concluded that one can not 
distinguish them from their decays, because they have the similar decay channels and the corresponding 
decay widths are comparable [12]. In this paper, we give the calculations of the decays + 
D'^j{2S) and S+ D'^ + Dj^(lD), we find that we can separate them, and according to the current 
experimental data of Belle Collaboration, we conclude that the new state Dfj{2700) is the first excited 
state Dfj{2S). These channels have also been considered in literatures, for example. Close et al [4, 23] 
give branching ratios and possible mixing of 2S and ID; Colangelo et al [7] assume the new state is the 
2S state, and extract the decay constant. 

The decay amplitude for — )• + D^j can be described in the naive factorization approach: 

T = ^V,sV,bai{D^p,\0){D''\J'^\B+), (1) 

where the CKM matrix element Vcs = 0.97334 and Vcb = 0.0415 [13], since we focus on the difference 
between D'^j{2S) and D~^j{lD), not on the careful study, we have chosen ai = ci + |c2 = 1, where ci and 
C2 are the Wilson coefficients [14] . We also delete other higher order contributions like the contributions 
from penguin operators. And 

{Dtj\JA^)=iF^M^+el, (2) 

Fv and are the decay constant and polarization vector of the meson Z)+ , respectively. 

We have solved the exact instantaneous Bethe-Salpeter equations [15] (or the Salpeter equations [16]) 
for 1~ states, the general form for the relativistic Salpeter wave function of vector 1~ state can be written 
as [17, 18]: 

+ APk{QA-) + - q± ■ el)f7{q±) + j^iPAii. - A± • 4)/8(g±), (3) 
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where the P, q and e^^ are the momentum, relative inner momentum and polarization vector of the vector 
meson, respectively; fi{q±) is the function of —q"^, and we have used the notation qj_ = q'^—{P-q/M'^)P'^ 
(which is (0, q) in the center of mass system). 

The 8 wave functions fi are not independent, there are only 4 independent wave functions [18], and 
we have the relations 

r r \ [q±f3iq±) + M'^f5{q±)] (mims - ujiu;2 + ql) ... f5{q±)M{-mim2 + ujiUJ2 + ql) 

M[mi + m2)qj_ (w-i — m2)qj_ 

f( \ [-9i/4(g±) + M2/g(g^)] (mia;2 - m2a;i) ^%{q^^M(mxU}2 - m2i>:{) 
J2[q±) = 777 \ .-"2 , J8{q±) = 7 r-2 • 

M(wi + uj2>qi [1^1 - ^2)qi 

In our method, the S — D mixing automatically exist in the wave function of 1~ state, because we 
give the whole wave function Eq. (3) which is = 1~, but some of the wave function are S wave, some 
are D wave, for example, see Figure 1-3, we show wave functions. One can sec that, for IS" and 25", the 
wave function /g and /e are dominate, and they are S wave, but there is little D wave mixing in these 
two states which come from the terms fs and /4. But for the third state, we labeled as ID state, the 
terms /a and /4 are dominate, but these two terms are not pure D wave, they will give contribute as a 
S wave [17]. So we conclude that the IS and 2S states are S wave dominate states, mixed with a little 
D wave (come from the terms /s and /4), while the ID state is a D wave dominate state (/s and /i), 
mixed with a valuable part of S wave (still come from the terms /s and /4). 

For cs vector 1~ state, our mass prediction for the first radial excited 2S state is 2673 MeV, and for 
ID, our result is 2718 MeV [18], so there are two states around 2700 MeV. 

In Ref. [18], we only give the leading order calculation for decay constant, which is Fy = '^■\fN~c J T^W/sl^) 



(2^) 

the whole equation should be 

and our results of the decay constants for vector cs system are: 

FvilS) = 353 ±21 MeV, 

Fv{2S) = 295 ±13 MeV, 

Fv{lD) = 57.1 ± 5.1 MeV. 

Where the uncertainties are given by varying all the input parameters simultaneously within ±5% of 
the central values in our model, and we will calculate all the uncertainties this way in this letter. The 
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Figure 1: wave functions of D*{1S). 
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Figure 2: wave functions of D*{2S). 
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Figure 3: wave functions of D*(1D). 

center values of S waves are little smaller than the estimates in Ref. [18], where ^^(15) = 375 it 24 MeV 
and Fv{2S) = 312 it 17 MeV, but the value of D wave decay constant is much smaller than the predict 
in Ref. [18] where we did not shown it, this is because we should not ignore the contribution from the 
term of /s when consider a D wave state. Our estimate of Fv{lS) = 353 it 23 MeV is close to the 
newer result Fd* = 268 ~ 290 MeV by Choi [19]. Our estimate of ^^(25") = 312 ± 17 MeV is a little 
larger than the estimate Fv{2S) = 243 it 41 MeV in Ref. [7], where they extracted it from the decay 
B+ ^ + D+,(2700) assuming the new state L>+j(2700) is the first radial excited state D^j{2S). 

According to the Mandelstam formalism [20], at the leading order, the transition amplitude for 
5+ — )• can be written as [21]: 

dq 



{D''\J^'\B^ 



(2vr)= 



DO rric + mu M B+ 



(5) 



where f is the recoil three dimensional momentum of the final state meson, ip^~^ is called the positive 
energy wave function, and (^^^^ = Toliy^ty^j^To- 

The wave function forms for pseudoscalar and are similar, for example, the wave function for 
B+ can be written as [22] 



+ : 



uju + ^b , i{mu - nib) , #70 i'^u + i^b) 

+ 7o 



' / i 7 ) (6) 

where a;^ = -\/m^"-i-g^ and a;;, = + g^; V'i(g), </'2(C) the radial part wave functions, and their 
numerical values can be obtained by solving the full Salpter equation of 0~ state [22] . 
The decay width is: 



where p^^^ and Mf^ are the three dimensional momentum and mass of the final new state D*'j ■ Mj^ 
come from the definition of decay constant in Eq.(2), but the square of polarization vector and transition 
amplitude of Eq.(5) are symbolized as X. 
So our result are: 

r{B+ D^D*+{2S)) = F^{2S) x (3.50 ± 0.38) x 10"^^ GeV, (8) 

T{B+ D'^D*+{1D)) = F^{ID) x (3.25 ± 0.30) x 10"^^ GeV. (9) 

If we ignore the mass difference between L'*+(25) and D*"'"(1D), and use 2700 MeV as input, the results 
become 

T{B+ D^D*+{2S)) = F^{2S) x (3.36 ± 0.25) x lO^^^ GeV, (10) 

V{B+ L>°D*+(1D)) = Fy(l-D) X (3.36 ± 0.25) x lO"^'' GeV. (11) 

So the difference between this two channel mainly come from the difference of decay constants. Then our 
predictions of branching ratios are: 

Br{B+ D^D*+{2S)) = (0.72 ± 0.12)%, (12) 

Br{B+ D^D*+{1D)) = (0.027 ± 0.007)%. (13) 

In Ref. [12], we have calculated the main decay channels of Dfj{2S) and D'^j{lD), and we have the 
following estimates: 

Br{D+j{2S) D^K+) = 0.20 ± 0.03 



and 

Br{D+j{lD) D^K+) = 0.32 ± 0.04, 

so we obtain: 

Br{B+ D°Djj{2S)) x Br{D+j{2S) D°K+) = (1.4 ± 0.5) x 10"^ (14) 

and 

Br{B+ D°Djj{lD)) x Br{D+j{lD) D^K+) = (0.9 ± 0.3) x 10"^ (15) 

Our estimate of Br{D^j{2S) D^K'^) ~ 0.20 is larger than than the estimate 0.11 in Ref. [4] and 
the estimate 0.05 in Ref. [5], if we use their value as input, we will obtain a smaller value of Br{B~^ — > 
D°D+j{2S)) x Br{Dfj{2S) D^K+). But our estimate of Br{D+j{lD) D^K+) ~ 0.32 consist with 
the value 0.34 in Ref. [5]. One can also see that, the final branching ratios depend strongly on the decay 
constants, but at this time few papers have calculated the values of decay constants for D^j(2S) and 
D'^j{lD). In Ref. [7], they assumed that the new state D^(2700) is D'^j{2S), and they extracted decay 
constant of D+j{2Sy. ^"^+^(25) = 243 ± 41 MeV. 
The Belle Collaboration have the data [3, 13] 

Br{B+ D°D+j{2700)) x Br{Djj{2700) D^K+) = (1.13 toU) x 10"^ 

Because we only calculated the decay widths of six main channels for Dfj{2S) {D'^j{lD)), and based 
on the summed width of these six channels, not full width, wc give a relative larger branching ratio 
Br{Dfj{2S) D^K~^) and Br{Dfj{lD) D^K~^), the real branching ratios should be smaller than 
our estimates, so our estimate of 5+ decay to D'^j{2S) is close to the data, while the estimate of 
decay to D^j{lD) is much smaller than the data, then we can draw a conclusion that the new state 
(2700) is Dtj{2S). 

We have another method to estimate the branching ratio, because the mass of is much heavier 
than the mass of D^, and the mass of Dfj{2700) is close to the mass of D*'^, as a rough estimate, we 
ignore the mass difference of Z'j}(2700) and Z'*+(2112), from Eq.(l) and Eq.(2), then we have 

Br{B+ ^ D^Dtj{2S)) ^ F^{2S) 

Br{B+ ^ DOD*s+{lS)) F^{1S)' ^ ' 

and from Particle Data Group [13], the branching ratio of B^ + Dfj{lS): 

Br{B+ D°D*+{1S)) = (7.8 ± 1.6) x 10"=^. 
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So we have: 



Br{B+ D'^D*+{2S)) ~ 



X (7.8 lb 1.6) X 10"^ ~ (5.4 ± 1.7) x 10 
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Br{B+ ^ D^D*+{1D)) 



Fum 

F^ilS) 



X (7.8 ± 1.6) X 10-3 



~ (2.0 ± 1.0) X 10 
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(18) 



This rough estimate results arc very close to our calculations (Eq.(12) and Eq.(13)), so we have the same 
conclusion that L»+j(2700) is D+j{2S). 

In Ref. [12], we estimate the full widths of Dfj{2S) and Dfj{lD) by six main decay channels, the 
estimated full widths are 46.4 ± 6.2 MeV for D+j{2S), 73.0 ± 10.4 MeV for Dfj{lD), comparing with 
experimental data, T = 112 ± 7 ± 36 MeV for £)+^(2690) and T = 108 ± 23^11 MeV for i:>+^(2700), there 
is still the possible that there are two states around 2700 MeV, one is the S wave dominate D^j{2S), the 
other is D wave dominate Dfj{lD). 

As summary, from the decays 5+ + D'^j{2S) and 5+ + Dfj{lD), we conclude that the 

new state D^(2700) is the first radial excited state D'^j{2S), and there may exist another state around 
2700, Djj{lD), with a mass around 2718 MeV, and a width 73.0 ± 10.4 MeV. 
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